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Abstract 

Metal nanostructures act as powerful optical antennasll|,|20 because collective modes of the electron 
fluid in the metal are excited when light strikes the surface of the nanostructure. These excitations, 
known as plasmons, can have evanescent electromagnetic fields that are orders of magnitude larger 



than the incident electromagnetic field. The 
tween plasmonically active nanostructuresy, 



argest field enhancements often occur in nanogaps be- 
4\, but it is extremely challenging to measure the fields 



in such gaps directly. These enhanced fields have app ications in surface-enhanced spectroscopies 
7], nonlinear opticslfl, ||3- ^ 



10], and nanophotonicsilllMlStl. Here we show that nonlinear tunnelling 



conduction between gold electrodes separated by a subnanometre gap leads to optical rectification, 
producing a DC photocurrent when the gap is illuminated. Comparing this photocurrent with low fre- 
quency conduction measurements, we determine the optical frequency voltage across the tunnelling 
region of the nanogap, and also the enhancement of the electric field in the tunnelling region, as a 
function of gap size. The measured field enhancements exceed 1000, consistent with estimates from 



surface-enhanced Raman measurements lll6l4l8h . Our results highlight the need for more realistic 
theoretical approaches that are able to model the electromagnetic response of metal nanostructures 
on scales ranging from the free space wavelength, A, down to ~ A/ 1000, and for experiments with 
new materials, different wavelengths, and different incident polarizations. 
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The probability of an electron being able to tunnel across a small gap decays exponentially 
with distance, so electron tunnelling can be used to measure distances with high precision, as 
happens in scanning tunnelling microscopy (STM). Moreover, nonlinearity in tunnelling conduc- 
tion can convert the optical frequency response of the metal into a DC current. This rectification 
process, originally proposed as a means of inferring the tunnelling time lll9n . allows the optical 
frequency voltage difference across the tunnelling gap to be determined. Evidence for microwave 
rectification ll20ll and optical rectification fc ll 12211 has been reported for STM experiments, but ther- 
mal and surface photo voltage effects [|20n complicate interpretation of the optical data. 

Suppose that the DC two-terminal current-voltage characteristics of some object are /(Vdc)> 
where I is the current and Vdc is the DC voltage bias. If the underlying conduction mechanism 
remains valid at some frequency scale 2co, then in the presence of some small AC voltage bias, 
Vac cos(cjt), in addition to Vbc, one expects Il20ll a total current of 



dl 1 d 2 I 

I » I(V DC ) + T^kc^Ac cos(wt) + -^|y DC ^A C cos 2 (wt) 
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In addition to /(Vdc)> m me presence of the AC excitation there is a rectified DC current 
proportional to the nonlinearity of the conductance and the square of the AC excitation amplitude. 
As demonstrated[20] using STM and microwave irradiation, one may use a low frequency AC 
voltage and lock-in amplifier to measure d 2 I/dV 2 using the second harmonic of the AC voltage. 
Using lock-in techniques to measure the DC current due to microwave irradiation of a known 
microwave voltage, Tu et al. Il20n find quantitative agreement between the rectified microwave 
current and the expectations of Eq. (OQ). 

In the optical case, in the limit that this rectification picture is valid, one would expect the 
radiation-induced AC response of the metal to lead to a photocurrent, J p hoto> in addition to the usual 
DC current, /(Vdc)- One may then infer the local, radiation-induced AC voltage, V opt , by com- 
paring /photo and detailed low frequency characterization of the conductance. In a tunnel junction, 
the differential conductance, dl/dV, at zero bias can be related to the interelectrode separation 
through the exponential decay of the conductance with increasing distance, G ~ exp(— f3(d — d )). 
Here d is the internuclear distance between closest atoms, do is the size of the metal lattice con- 
stant, and (3 is the attenuation factor. The average local electric field at the closest interelectrode 
separation, V opt /(d — d ), may then be inferred and compared with the incident field known 
from the incident optical intensity, giving a quantitative measure of the field enhancement due to 



the metal optical response. The V opt inferred experimentally is that experienced by the optical- 
frequency tunnelling electrons. Since both V opt and d — d are found from tunnelling, a very local 
process, this is self-consistent. Voltage drops within the metal far from the tunnelling gap are not 
relevant. However, there are subtleties to consider, including the role of screening in the metal and 
the locality of the tunnelling process. We discuss this further in Supplemental Information. 

A quantum treatment of radiation interacting with a nanoscale junction involves photon-assisted 



tunnelling Q23Q and has been applied to atomic-scale metal contacts 11241 12511 . When a junction 
is illuminated with radiation of energy hui, the junction's plasmonic response results in an AC 
voltage, V opi at frequency u across the junction. Following the Tien-Gordon approach [126ll . the 



quantum correction to the DC current induced by the radiation in the limit of small AC amplitudes 

(eKpt < tiw) is given by 



/(Vbc> Vopt, of) — /(Vbc) — ^Kpt 



/(Vbc + Tiu/e) - 2I(V DC ) + I(V DC - tkj/e) 



(2) 



(hw/e) 2 

If the tunnelling nonlinearity is small on the voltage scale fko/e, this expression reduces to Eq. © 
with V op t playing the role of Vac- This occurs when the transmittance of the junction r(E) is 
smooth as a function of energy near the Fermi level of the leads, E F , throughout E F ± ftw. 

Details of the sample fabrication and the optical and electronic measurements are described in 
Methods. Nanogaps are formed in the Au constrictions via electromi gration H27H . The electromi- 
gration process is stopped once the zero bias DC conductance of the gaps is less than 2e 2 / h = G . 
Typical conductances in this experiment are on the order of 0.6 Go to 0.016 Go indicating average 
gaps, d — do, ranging from 0.03 nm to 0.23 nm. Figure 1 shows the measurement system and basic 
characterization of a representative nanogap. The illumination wavelength is 785 nm, with a peak 
intensity of 22.6 kW/cm 2 , and measurements are performed in vacuum at 80 K. 

The DC I — V characteristics of the samples after electromigration are weakly nonlinear, with 
no sharp features, as expected for clean vacuum tunnel junctions. These junctions, ideally devoid 



of molecules, show no molecular surface-enhanced Raman scattering signal 1116111711 . as expected. 
Fig. 2 shows representative d 2 I/dV 2 vs. Vbc curves and simultaneously recorded I p hoto vs - Vbc 
data for three different samples, acquired with the beam fixed over the nanogap. For each set of 
curves shown, the amplitude of Vac nas been set such that jV£ c d 2 I/dV 2 and J p hoto are the same 
amplitude. The fact that J p hoto is directly proportional to d 2 I /dV 2 consistently over the whole DC 
bias range, including through the Vbc where both change sign, strongly implies that the photocur- 
rent originates from the rectification mechanism. In this case, when \Vl c d 2 I /dV 2 = /photo, this 



implies Vac = Vopt, where V opt is the amplitude of the optically induced voltage between the two 
electrodes that oscillates at uj/2% = 3.8 x 10 14 Hz. When the rectification mechanism is responsi- 
ble for the photocurrent, comparison with low frequency electrical measurements gives a means of 
quantitatively assessing the plasmon-induced AC voltage across the nanoscale interelectrode gap. 
As a consistency check, we have mapped the optically rectified current, J p hoto> at fixed Vbc (while 
monitoring /, dl/dV, and d 2 I /dV 2 , to ensure junction stability during the mapping procedure). 
As shown in Fig. le, J p hoto is only produced when the gap itself is illuminated. 

Previous tunnelling optical rectification experiments have been hampered by thermal voltages 
and thermal expansion. We do not expect thermal voltages in our experiment since the device 
is a comparatively symmetric design, and with a centered laser spot there should be no tempera- 
ture gradient between the electrodes, even in the presence of minor structural asymmetries. The 
macroscale unilluminated parts of the electrodes also act as thermal sinks. Differential thermal 
voltages caused by asymmetry of illumination would show up as systematic trends in /photo as a 
function of laser position, and Fig. le shows no such trends. Thermal expansion should be irrel- 
evant, since the electrodes are intimately mechanically coupled to the substrate, with essentially 
no unsupported electrode material. To check for thermal effects, we looked for changes in dl/dV 
at uichopper, since any laser-induced expansion should modulate the interelectrode gap and thus the 
conductance. No measurable effect was found to better than the 10~ 4 level. 

Optical rectification is further supported by measurements of J p hoto vs. incident laser power. In 
Fig. 3 we see the expected linear dependence. The best fit line has a slope of 0.37 nA/(kW/cm 2 ). 
Further, as expected, Vbc can be selected such that d 2 I/dV 2 is zero (not always at Vbc = 0), and 
at that bias J p hoto measurements show no detectable current. 

The validity of the rectification picture is not completely obvious for this system a priori. De- 
tailed calculations of the transmission as a function of energy, as shown in Fig. 4, demonstrate 
that for this particular material and illuminating wavelength, corrections to the classical rectifica- 
tion picture should be relatively small. At higher DC biases (> 100 mV, not shown), we do often 
observe deviations of the photocurrent from the measured [1 / A)Vl c d 2 1 / dV 2 . 

To infer the nanogap size from the conductance measurements, we have performed theoretical 
simulations of the breaking process of gold atomic contacts. In these simulations, we have used 
density functional theory to determine the geometry of the junctions as well as to compute the 



transport properties within the Landauer approach, as detailed in Ref. Il28h . In this approach the 



low-temperature linear conductance is given by G = G r(E F ), where r(E F ) is the transmission 



of the contact evaluated at the Fermi energy, E-p. To simulate the rupture of the gold nanowires, 
we start out with a single-atom contact {e.g., Fig. 4a). We separate the electrodes stepwise, re- 
optimize the geometry in every step, and compute the corresponding conductance. As in Fig. 4c, 
in the relevant range, the linear conductance decays exponentially with an attenuation factor of 
/3 ~ 1.8 A -1 . Extensive simulations (Supplementary Information) show that this exponent varies 
within 10% depending on the geometry of the contact. 

Using — 1.85 A -1 (averaged over different junction structures), the interelectrode gap, d—d , 
can be found from the measured dl/dV. With this and the inferred V opt , the enhanced electric 
field local to the tunnelling gap may be determined. The experimental uncertainty in V opt is ap- 
proximately 5%, based on the comparisons shown in Fig. 2. The experimental uncertainty in 
the tunnelling conductance at zero bias is ~ 1%. The dominant uncertainty in the inferred inter- 
electrode separation is systematic, due to variation of (3 with morphology of the metal surfaces. 
Assuming that d — d is the appropriate scale for computing the local electric field, for the samples 
in Fig. 2, the optical electric fields (RMS) are 2.1 x 10 s V/m, 5.7 x 10 8 V/m, and 3.6 x 10 8 V/m, re- 
spectively. Based on gaussian beam shape and incident power, the free-space incident RMS field 
is 2.9 x 10 5 V/m. Therefore we determine local enhancement factors of the electric field to be 
718, 1940, and 1230. These very large enhancements exceed expectations from idealized FDTD 



calculations II 1611 for wider gaps , and are consistent with observations of single-molecule Raman 



emission in such junctions 
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18Q. At sufficiently high electric field strengths the gap geometry 
can become unstable. We believe that this is the reason why I — V curves of some devices exhibit 
instabilities. In the absence of the optical field or at considerably reduced incident laser fluences, 
the I — V curves are very stable. 

Our calculations clarify why the classical rectification formula works in our case. The transmis- 
sion function (Fig. 4d) of these gold contacts shows only a weak energy dependence around the 
Fermi energy in the interval [E F — ftw, E F + ftw). This indicates that the quantum result of Eq. © 
reduces to the classical rectification formula of Eq. ©. This smoothness is because transport in 
gold contacts is dominated by the s atomic orbitals, with a relatively flat density of states. Similar 
physics should therefore be expected in other s-valent metals, such as Ag. 

In several devices, conductance decreased over the course of hours of measurements, as the 
metal electrodes rearranged themselves, allowing us to track variations in V opt (and inferred en- 
hanced electric field) as a function of interelectrode separation (Fig. 5). Despite device-to-device 
variability, the general trend shows increasing enhancement with decreasing d — d . This trend is 



slower than l/(d — d ) and weaker than classical expectations Il29ll . Such a comparatively weak 



increase in enhancement at low d — do is qualitatively consistent with calculations that examine 



30D. How- 



quantum tunnelling corrections to plasmonic enhancements in nanoparticle dimersH29, 
ever, those calculations predict an actual decrease in enhancement when d — d is small enough 
that interelectrode conductance is strong enough to short out the nanogap plasmons. From our 



data and previous observations Hl7H . this decrease in enhancement does not happen in gold nano- 
j unctions until interelectrode conductance exceeds Go- Our measurements highlight the need for 
detailed, realistic theoretical treatments of these nanostructures, incorporating optical interactions, 
quantum effects, and dynamical screening. 
Methods 

Samples are fabricated using standard electron-beam lithography techniques to define metal 
(1 nm Ti/15 nm Au) constrictions 100 nm in width and 600 nm in length between larger contact 
pads, on highly doped n-type Si substrates covered with 200 nm of thermally grown Si02- After 
cleaning in oxygen plasma, devices are wirebonded to a chip carrier, installed in a microscope 
cryostat, and cooled to 80 K. 

After electromigration the sample is illuminated with a 785 nm laser focused to a gaussian spot 
with FWHM 1.9 /im and a peak intensity of 22.6 kW/cm 2 . The laser is rastered over the sample 
surface and the Raman response of Si substrate is measured at 520 cm -1 . The Au electrodes 
attenuate the Si Raman emission, allowing precise mapping of the nanogap structure as shown in 
Fig. Id. Once the nanogap is located via the Si Raman map, the laser is centered over the gap. 
The laser polarization is aligned along the original constriction. 

The sample is electrically characterized via the apparatus presented in Fig. la using two 
lock-in amplifiers and a current amplifier. The first lock-in amplifier applies an AC excitation 
(Vac cos wit) to the sample at uj\ = 2.0 kHz and simultaneously measures the first (oc dl/dV) 
and second harmonic response (= \V^ c d 2 1 / dV 2 ) of the nanogap to this excitation. There is no 
measurable change in either of these responses when the laser illumination is added or removed. 
The first harmonic response is essentially in phase with the oj\ reference. The second harmonic 
response is phased relative to the 2u>i reference by 180 degrees as defined by the lock-in, such 
that the integrated d 2 I/dV 2 signal gives back the dl/dV signal. The second lock-in amplifier is 
referenced to (and in phase with) an optical chopper at 0J c h oppe r = 232 Hz and measures the current 
generated due to the radiation (7 p hoto)- A summing amplifier is used to add a DC voltage Vdc to 
Vac- With this arrangement, 7 (Vdc), dl/dV, d 2 I/dV 2 , and 7 photo as a function of Vdc can be 
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measured simultaneously. 

[Supplementary Information] is linked to the online version of the paper at www.nature.com/nature. 
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FIG. 1. Measurement approach and layout, a) Schematic of electrical characterization measurement, b) 
Colorized SEM image of typical nanogap device. The actual tunnelling gap is not resolvable by SEM. 
c) Larger SEM view of electrodes and nanogap. d, e) False-colour maps of Si Raman line (d) and the 
photocurrent (e) for the device shown in (c). In (d) blue indicates a weak Si signal, thus indicating the 
location of the Au electrodes. In (e) the color bar runs from red (20 nA) to blue (0 nA). The photocurrent is 
clearly localized to the nanogap, to within optical resolution. 
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FIG. 2. Demonstration of optical rectification, a-c) Photocurrent (red, /photo) an d \Vj* c d 2 I/dV 2 (black) as 
a function of Vbc f° r three different samples. The shared shapes of these curves, including changes of sign, 
demonstrate that the photocurrent arises from the rectification process. Note the differing current scales for 
the three devices. The different conductances and nonlinearities presumably result from microscopic details 
of the electrode surfaces and geometry. Insets: Conductance in units of /xA/V for each device. Gap distances 
from tunnelling measurements are 1.4 A, 0.44 A and 0.92 A. Inferred V opt values are 30 mV, 25 mV, and 
33 mV, with uncertainties of 10%. Inferred electromagnetic fields are 2.1 x 10 8 V/m, 5.7 x 10 8 V/m and 
3.6 x 10 8 V/m yield field enhancements of 718, 1940, and 1230, respectively. 
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FIG. 3. Further evidence for optical rectification. The rectified photocurrent, 7 p hoto> as a function of the 
incident laser power. The error bars indicate one standard deviation of the I p hoto measurements at each 
laser power. The linear power dependence (see text) was found by a weighted least squares fit (red line). 
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FIG. 4. Theoretical basis for validity of rectification, a) Geometry considered in the theoretical analysis 
of the distance dependence of the linear conductance. We start with a single-atom contact grown along 
the (111) direction, b) Stretched geometry with a distance d between the gold tips, c) Calculated linear 
conductance as a function of the interelectrode distance (full circles). The solid line shows the fit to the 
exponential function G = Go exp[— /3(d — do)]. The fit parameters are indicated in the graph, d) Zero-bias 
transmission as a function of the energy for the different geometries of panel (c), with d increasing from 
top to bottom. Note the logarithmic vertical axis, and that the transmission remains a smooth function of 
energy over a broad range around the Fermi level. 
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FIG. 5. Field enhancement (left axis) at the tunnelling region as a function of gap distance (top axis) 
for five devices (shown in different colours) measured a number of times at 80 K; the gap distance (and 
therefore the enhancement) changes as the gap geometry evolves over the course of repeated measurments. 
The enhancement and gap distance are derived from measurements of the conductance (bottom axis shows 
ln(Go/G)) and the optical frequency voltage (V^ pt , right axis) based on the assumption that d — do is the 
relevant scale over which V^pt falls, and /3 = 1.85 A -1 . Error bars indicate the statistical uncertainty in 
V opt and In G. 
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Supporting Information: 
Optical Rectification and Field Enhancement in a Plasmonic Nanogap 

I. TUNNELLING AND FIELD DISTRIBUTIONS IN NANOGAPS 

The experimentally extracted parameter in the rectification measurement is the optical fre- 
quency potential difference, V opt , experienced by the tunnelling electrons. To go from this poten- 
tial difference to an estimate of the enhanced electric field within the tunnelling gap, it is necessary 
to estimate the distance over which this voltage is dropped. In the main manuscript we use d — d , 
the interelectrode spacing inferred from the (low frequency) tunnelling conductance. The question 
arises, is this reasonable given that V opt is oscillating at optical frequencies? There are two relevant 
issues: the local character of the tunnelling process, and the penetration of electric field (and hence 
voltage drops) into the metal electrodes. 

The rectification process is sensitive to the chemical potential difference experienced by pre- 
cisely the tunnelling electrons. Because of the local nature of the tunnelling process, the chemical 
potential difference we infer experimentally is measured between the "initial" and "final" positions 
of the tunnelling electrons. The inferred V opt is the voltage difference that drives the (extremely 
local) tunnelling process. It is self-consistent to use the gap size extracted from tunnelling con- 
ductance as the relevant distance scale, since the local character of the tunnelling electrons is not 
expected to be frequency dependent (that is, there is no expectation that tunnelling takes place over 
longer distances at higher frequencies). Additional voltage drops on longer scales within the metal 
electrodes are perfectly allowed, but not relevant to the inferred V opt since the local tunnelling 
process is not sensitive to them. The observed consistency between J p hoto and d 2 I/dV 2 over the 
whole bias range would be very difficult to understand if the local nature of the tunnelling pro- 
cess was fundamentally different between low and optical frequencies. This is certainly a subtle 
issue. The definition of local electrostatic potential and its relationship with chemical potential 
in nanoscale junctions is not simple si , since electronic distributions are driven out of equilibrium 
on a distance scale comparable to the inelastic scattering length. A full treatment of optical- 
frequency tunnelling with realistic electronic structure of the contacts (well beyond the scope of 
the present work) would be needed to resolve this issue unambiguously. The experimentally deter- 
mined quantities in the main manuscript (tunnelling conductance, V opt ) could be compared with 
more sophisticated calculations in the future. 
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Regarding field penetration into the metal, care is also required. The skin depth, which de- 
scribes the characteristic attenuation length for the electric field of a plane wave normally incident 
on a metal surface, is related to the complex dielectric function at the incident frequency. For gold 
at the incident wavelength of 785 nm, the skin depth is approximately 25 nm. However, the skin 
depth calculation ignores collective electronic response such as plasmons. In fact, in plasmonic 
nanostructures, while the field due to the incident radiation does penetrate the metal on length 
scales comparable to the skin depth, the total electric field including that produced by the plasmon 
response of the electron fluid can be confined much more strongly than this. 

Supplemental Figures ISTHS31 demonstrate this explicitly. These figures show the results of a 
finite element simulation (using the commercial package COMSOL) of a model nanogap structure 
illuminated at the nanogap plasmon resonance. The figures show top views of the metal structure 
at three different magnifications indicated by the scale bars. The metal is modeled as Au (using 
the commonly employed Johnson and Christy dielectric function 32 ) 25 nm in thickness (chosen 
for computational efficiency). The nanogap separation is 10 nm, except in the center of the device, 
where a protrusion on the left-hand electrode reduces the gap to 5 nm. No dielectric substrate 
is included, but other models including such a substrate give qualitatively identical results. The 
incident radiation is assumed to be a plane wave at normal incidence with a free space wavelength 
of 850 nm and a polarization directed along the line connecting the two electrodes (left-right in the 
figures). Upper figures show the magnitude of the total electric field (normalized by the incident 
amplitude) on a linear color scale, while lower figures use a (natural) logarithmic color scale. 

It is immediately clear that the vast majority of the potential drop between the two electrodes 
takes place across the vacuum at the center of the nanogap. Electric field does extend into the 
metal, but it is smaller than the nanogap field, by a factor of more than 50 in this case, a larger 
difference than one would expect based on a skin depth analysis that ignores plasmonic effects. 
For smaller electrode separations, this disparity grows even further, with essentially 100% of the 
(plasmon-dominated) optically induced potential difference dropping across the vacuum. 

The situation is, of course, complicated in the real devices by the true electronic structure 
of the metal surfaces at the atomic scale, and tunnelling conduction. For very small gaps and 
distance scales, one must not trust the quantitative accuracy of finite-element or finite-difference 
time-domain calculations based on idealizations (infinitely sharp dielectric interfaces; bulk- like 
dielectric functions for the electrons and background charges at interfaces). 

However, in plasmon resonant structures such as these, the fact remains that a very large frac- 
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tion of the potential drop takes place between the electrodes, and the skin depth is not the relevant 
length scale to describe electric field confinement. The excitation of plasmon collective modes 
strongly perturbs the charge density extremely local to the surface, leading to field distributions 
that can differ greatly from those inferred in the absence of such collective effects. For example, 
Liebsch S3,S4 shows that the presence of a surface plasmon excitation in a silver film or nanopar- 
ticle dynamically leads to surface charge density within a couple of atomic radii of the surface, 
with a length scale determined by quantum mechanical electronic structure effects and not by 
macroscopic electrodynamic effects such as the classical skin depth. 

II. THEORETICAL METHODS 

Here we describe in more detail the theoretical methods used for the computation of the struc- 
tural, electronic, and electrical properties of gold nanocontacts. Moreover, we present an analysis 
of the geometry dependence of the conductance decay with interelectrode distance. 

A. Electronic structure and contact geometries 

For the determination of the electronic structure and geometry of the gold contacts we employ 
density functional theory (DFT) as implemented in the RI-DFT module of the quantum chem- 
istry package Turbomole 5.9 S5 > S6 . In all our calculations we use the open-shell formalism, the 
exchange-correlation functional BP86 S7 ' S8 , and Turbomole's standard Gaussian basis set "def- 
SVP" of split- valence quality with polarization functions S6 > S9 > S10 . Total energies are converged 
to a precision of better than 10~ 6 atomic units, and structure optimization is carried out until the 
maximum norm of the Cartesian gradients has fallen below values of 10~ 4 atomic units. 

An example of the contact geometries considered in this work is shown in Fig. IS4K. which 
corresponds to the junction of Fig. 4 in the manuscript. In this case, a single-atom contact is 
modeled by means of the extended central cluster (ECC) with 63 atoms displayed in Fig. IS4T3. It 
consists of two pyramids grown along the (111) direction of a fee lattice, which contain, beside 
a common tip atom, four layers with 3, 6, 10, and 12 atoms. We formally divide the contact 
geometry into left (L), central (C), and right (R) parts, and relax the atomic positions of the atoms 
in the C region, i.e. those of the tip atom and the neighboring layers with three atoms. In contrast, 
the atoms in the L and R regions are kept fixed, and the lattice constant is set to 4.08 A, the 
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experimental bulk value for gold. 

In order to simulate what happens when the wires are broken and when one enters the tunnelling 
regime, we stretch the ECC (Fig. |S4B) by separating the L and R regions in steps of 0.4 a.u. 

0.21 A). After every stretching step, the geometries are reoptimized by relaxing the C region. 



B. Conductance calculations 



To compute the charge transport properties we apply a method based on standard Green's func- 
tion techniques and the Landauer formula expressed in a local nonorthogonal basis, as described 
in Ref. Q28[|. Briefly, the local basis allows us to partition the basis states into L, C, and R ones, ac- 
cording to the division of the contact geometry (Fig. |S4l . Thus, the Hamiltonian (or single-particle 
Fock) matrix H, and analogously the overlap matrix S, can be written in the block form 



H 



/ 



V 



Hll H L c 
Hc*L Hcc Hcr 
Hrc Hrr 



\ 



(3) 



Within the Landauer approach, the low-temperature conductance is given by G = Got(Ef), where 
Go = 2e 2 /h and r(Ep) is the transmission of the contact evaluated at the Fermi energy. The 
energy-dependent transmission r(E) can be expressed in terms of the Green's functions as si : 



t{E) = Tr [T L G r cc T R G 



cc\ > 



(4) 



where the retarded Green's function is given by 



G r cc (E) = [EScc - H cc - YT L {E) - ^(E)]- 1 



and G cc = [G cc ]^. The self-energies adopt the form 



(5) 



^x(E) = (Hex — ES C x) gxx(E) (Hxc — ES X c) 



(6) 



the scattering rate matrices are given by T X (E) = — 21m [S^- (£?)], and g r xx (E) = (ESxx ~ 
Hxx) -1 are the electrode Green's functions with X = L, R. 

In order to describe the transport through the contact shown in Fig. IS4R. we first extract Hc^ 
and Sec an d the matrices Hex an d Sex from a DFT calculation of the ECC in Fig. IS4B. The 
blue-shaded atoms in regions L and R of Fig. IS4R are assumed to be those coupled to the C 
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region. Thus, Hex and Sex, obtained from the ECC, serve as the couplings to the electrodes in 
the construction of H r x (E). 

On the other hand, the electrode Green's functions g r xx (E) in Eq. © are modeled as surface 
Green's functions of ideal semi-infinite crystals. In order to obtain them, we have computed 
separately the electronic structure of a spherical fee gold cluster with 429 atoms. By extracting the 
Hamiltonian and overlap matrix elements connecting the atom in the origin of the cluster with all 
its neighbors and by using these as "bulk parameters", we construct a semi-infinite crystal which 
is infinitely extended perpendicular to the transport direction. The surface Green's functions are 
then calculated from this crystal with the help of the decimation technique S11 > S12 . In this way we 
describe the whole system consistently within DFT, using the same nonorthogonal basis set and 
exchange-correlation functional everywhere. 

We assume the Fermi energy Ep to be fixed by the gold leads. From the AU429 cluster we obtain 
E F = —5.0 eV. Let us point out that the Fermi energy is used only to read off the conductance 
from the transmission function r(E), but it does not enter into our calculations otherwise. 

C. Geometry dependence of the conductance decay with interelectrode distance 

When an atomic contact is broken, the conductance depends exponentially on the interelectrode 
distance d, i.e. G = G exp(—/3(d — d )). Here, d describes the prefactor of the exponential law 
and (3 is the attenuation factor. Both constants depend on the surface morphology, in particular on 
the local atomic environment at the interface region as well as on the crystallographic direction 
along which the contact is grown. We have studied the influence of these two factors by performing 
different DFT-based simulations of the breaking process of gold nanowires. 

In Fig. IS5A-D we show the results obtained for a contact grown along the (111) direction. 
As compared to the geometry considered in Fig. 4 of the manuscript, the starting geometry is a 
one- atom-thick contact with a dimer (or two-atom chain) in the narrowest part. We have stretched 
this contact following the procedure described above (see paragraph "Electronic structure and con- 
tact geometries"), and we have computed both the linear conductance and the energy-dependent 
transmission at the different stages of the elongation process. As it is shown in Fig. |S5C, in 
the range relevant for the experiments, the conductance decays exponentially with an exponent 
= 1.93 A" 1 . This value is approximately 7% larger than the value obtained for the contact of 
Fig. 4. Furthermore, in Fig. IS5t > one can see that the transmission is still a smooth function of the 
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energy in a range larger than 1 eV around the Fermi energy. 

In Fig.lS5E-H we present the corresponding results for a similar contact with a dimer structure, 
but this time grown along the (100) direction of the fee lattice. In this case (3 takes a value of 
1.88 A -1 , which is 4% larger than those found in Fig. 4. Again, the transmission is a rather 
featureless function of the energy in the relevant conductance range. 

To conclude, we have analyzed different types of gold atomic-sized contacts. Our results indi- 
cate that the attenuation factor (3 varies roughly within 10% of the value reported in the manuscript. 
Additionally, they corroborate that the transmission is a rather smooth function of the energy in an 
interval of approximately 2 eV around the Fermi energy for the tunnelling regime explored in the 
experiments. 
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FIG. S 1 . Finite element simulation of the field distribution in a gold nanogap structure illuminated at its 
nanogap plasmon resonance. Model parameters are described in the Supplemental Information text. A 
spatial scale bar is shown. The upper figure shows the magnitude of the local electric field normalized by 
the magnitude of the incident field on a linear color scale. Maximum field enhancement is by a factor of 
100, at the closest interelectrode separation. The lower figure shows the same data plotted with a (natural) 
logarithmic color scale. 
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FIG. S2. A closer view of the model of Fig.[ST] While there is some electric field penetration into the metal, 
the vast majority of the total electric field is confined to the vacuum region between the electrodes. That 
is, almost none of the potential difference between the electrodes is dropped within the metal. Note further 
that the tunnelling process is sensitive only to the potential difference extremely local to the tip atoms. 
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FIG. S3. A closer view of the model of Fig.[S2] Again, while there is some electric field penetration into the 
metal, the vast majority of the total electric field is confined to the vacuum region between the electrodes. 
That is, almost none of the potential difference between the electrodes is dropped within the metal. The 
skin depth parameter is not relevant in determining this, because the skin depth does not account for fields 
produced by collective plasmon response of the electron liquid. 
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FIG. S4. Quantum transport scheme. The atomic-sized contact (A) is divided into a C region and two 
semi-infinite L and R electrodes. Using a similar division for the ECC (B), information on the electronic 
structure of the C region (Hccs Sec) as well as the CL and CR couplings (Hcl, &cl and Hcr, &cr) is 
extracted. The electrode surface Green's functions g r LL and g T RR , needed for the computation of the self- 
energies "E r L and £ j, [Eq. ©], are determined in a separate calculation as explained briefly in the text and 
in detail in Ref. Q28D- 
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FIG. S5. (A) Starting geometry used to study the distance dependence of the conductance of a gold contact. 
This geometry corresponds to a one-atom-thick contact grown along the (111) direction of a fee lattice with 
a dimer structure (or chain of two atoms) in the middle. (B) Stretched geometry with the distance d between 
the gold tips. (C) Calculated linear conductance as a function of the interelectrode distance (full circles). 
The solid line shows the fit to the exponential function G = Gq exp{—(3(d — do)]. The tit parameters do and 
(3 are indicated in the graph. (D) Zero-bias transmission t(E) as a function of the energy (measured with 
respect to the Fermi energy) for the geometries of panel C. (E-H) The same as in panels A-D for a dimer 
contact grown along the (100) direction. 
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